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A light front formalism for deep inelastic lepton scattering from finite nuclei is developed. In 
particular, the nucleon plus momentum distribution and a finite system analog of the Hugenholtz- 
van Hove theorem are presented. Using a relativistic mean field model, numerical results for the plus 
momentum distribution and ratio of bound to free nucleon structure functions for Oxygen, Calcium 
and Lead are given. We show that we can incorporate light front physics with excellent accuracy 
while using easily computed equal time wavefunctions. Assuming nucleon structure is not modified 
in-medium we find that the calculations are not consistent with the binding effect apparent in the 
data not only in the magnitude of the effect, but in the dependence on the number of nucleons. 



I. INTRODUCTION 



The nuclear structure function F2A ( x ) is smaller than A times the free nucleon structure function AF2N {x) for values 
of x in the regime where valence quarks are dominant. This phenomenon, known as the European Muon Collaboration 
(EMC) effect Jl[], has been known for almost twenty years. Nevertheless, the significance of this observation remains 
unresolved even though there is a clear interpretation within the parton model: a valence quark in a bound nucleon 
carries less momentum than a valence quark in a free one. There are many possible explanations, but no universally 
accepted one. The underlying mechanism responsible for the transfer of momentum within the constituents of the 
nucleus has not yet been specified. One popular mechanism involves ordinary nuclear binding which, in its simplest 
inculcation, is represented by evaluating the free nucleon structure function at a value of x increased by a factor of 
the average separation energy divided by the nucleon mass e/M ~ 0.04. The validity of this binding effect has been 
questioned; see the reviews |g 0, ||. 

The Bjorken variable x is a fraction of the plus component of momentum, and the desire to obtain a more precise 
evaluation and understanding of the binding effect lead us to attempt to obtain a nuclear wave function in which 
the momentum of the nucleons is expressed in terms of this same plus component. Therefore we applied light front 
dynamics to determining nuclear wave functions || . In this formalism one defines x ± = x° ± a; 3 and quantizes on 
equal x + surfaces which have a constant light front time, r. The conjugate operator P~ acts as an evolution operator 
for t. The plus momentum is canonically conjugate to the spatial x~ variable. This light front formalism has a variety 
of advantages J?], ||, ||, [l(], [nl 12. ll| 14, |l5) and also entails complications [[l6| 



Our most recent result ]fl7f is that the use of the relativistic mean field approximation, and the assumption that 
the structure of the nucleon is not modified by effects of the medium, to describe infinite nuclear matter leads to no 
appreciable binding effect. The failure was encapsulated in terms of the Hugenholtz-van Hove theorem [^J which 
states that the average nuclear binding energy per nucleon is equal to the binding energy of a nucleon at the top 
of the Fermi sea. The light front version of this theorem is obtained from the requirement that, in the nuclear rest 
frame, the expectation values of the total plus and minus momentum are equal. The original version of the theorem 
was obtained in a non-relativistic theory in which nucleons are the only degrees of freedom. Here, the mesons are 
important and the theory is relativistic, but the theorem still holds. This theorem can be shown to restrict O] the 
plus momentum carried by nucleons to be the mass of the nucleus, which in turn implies that the probability for a 
nucleon to have a plus momentum k + is narrowly peaked about k + — Ma/A = M. Thus the only binding effect 
arises from the average binding energy, which is much smaller than the average separation energy. Therefore dynamics 
beyond the relativistic mean field approximation must be invoked to explain the EMC effect. This conclusion was 
limited to the case of infinite nuclear matter, and the computed nuclear structure function could only be compared 
with data on finite nuclear targets extrapolated to the limit A — > 00. The goal of the present work is to extend the 
results to finite nuclei; the main complication arises from the spatial dependence of the nucleon and meson fields. 

We briefly outline our procedure. In Sections II and HI we present the covariant deep inelastic scattering formalism 
of Ref. [^0| and derive its representation in terms of nucleon single particle wave functions. The plus momentum 
distribution follows from this representation in Section IV where we also derive new version of the Hugcnholtz-van 
Hove theorem. Then we present the results of analytic and numerical calculations in Section the latter giving an 
A dependence of the ratio function contrary to experimental results. This again gives the result that the use of the 
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FIG. 1: Feynman diagram for deep inelastic scattering. A nucleus of momentum P is struck by a virtual photon of momentum 
q. We label nucleon momentum k, and quark momentum p. 



relativistic mean field approximation, combined with the assumption that the nuclear medium does not modify the 
structure of the nucleon, cannot describe the EMC effect. The reasons for the subtle differences between the results for 
finite nuclei and nuclear matter are detailed in Section VI. Finally, we summarize and discuss possible implications. 



II. NUCLEON GREEN'S FUNCTION FOR FINITE NUCLEI 



We begin with the covariant plus momentum distribution function 



Mv) 



d 4 k 



k° 



M 



Tr 
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2P+A 



X A {k,P) 



where we identify 

X A {k 1 P) = -if d 4 x J d 4 ye~ lk - (x - y) G c (x,y), 

where G (x,y) is the connected part of the nucleon Green's function: 

iG(x,y) = (P\T+{i;'(x)i?(y)}\P). 



(2.1) 



(2.2) 



(2.3) 



This result is directly determined from the Feynman diagram in Fig. [l] following Ref . [^0) . So far this is independent of 
the particular relativistic mean field model, but for concreteness we use a Quantum Hadrodynamics (QHD) Lagrangian 
p2| , p3[ , specifically QHD-I as in Ref. pl[ , where the nucleon fields, ip', that appear in Eq. ( |2.3| ) are those appearing 
in the Lagrangian. Light front quantization requires that the plus component of all vector potential fields vanishes, 
and this is obtained by using the Soper-Yan transformation pi, p5| 



^'{x) = e" l9 " A(a;) V(x), 



d + A{x) = V + {x). 



(2.4) 



to define the nucleon field operator ip for various models [ p~7[ . This transformation allows the use of the eigenmode 
expansion for the ip fields which have been obtained previously in Ref. |2l[] 



ip(x) = [a a u a {x) + b\, v a (x)] 

a 



(2.5) 



where a a and are (anti-)nucleon annihilation operators and we define z = —x~/2 with d + = 29_ = ~d z and 
(x±,z) = x which allows us to treat the minus and perpendicular coordinates on equal footing. The u a and v a are 
coordinate space 4-component spinor solutions to the light front Dirac equation with eigenvalues p^/2 = M — e a . 
To simplify the analysis we will temporarily ignore electromagnetic effects, but we will include them in the final 
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numerical results. The light front mode equations in QHD-I are obtained by minimizing the P operator (light front 
Hamiltonian) with the constraint that P + = P~ . The result is 



-idz\i> a ) 



[a± ■ (P± - 9vV ± ) + 0(M + g s cf>)] |^+) 
[-id, + 2g v V-] \1>+) 

+ [ a± ■ (p± - g v V ± ) + f3(M + g s 4>)} \^~), 



with 



u a (x) = V2P+A (x\i> a ) 

a±|^) = ^yVliW = h£> 

d+V^ = d + V^-d t "V + . 



(2.6) 
(2.7) 

(2.8) 
(2.9) 
(2-10) 



Using standard manipulations |23| and defining sp as the energy of the highest occupied state, we find the Green's 
function to be 



G(x,y) = J2 u «(x)u a (y)e- l9AA{x) - Aiv)] I ^ e -^t* + -« + )/2 



dk- 
~2x 



k - p a + ie 



— + 2m5(k -p a )e(e F -e a ) 



= G D (x ) y)+G u (x,y) ) 



(2.11) 



where the superscripts D and C represent the disconnected and connected parts of the nucleon Green's function, 
respectively. The connected part is relevant to deep inelastic scattering and is given by 



G c (x 1 y) = i 2J u a {x)u a (y)t 



-igvlA{x)-A(y)] -ip-(x+-y+)/2 



(2.12) 



where the sum is over occupied levels a in the Fermi sea F. We now substitute Eq. (2.12) into Eq. (2.2), first defining 
(fej_, k + ) = k where k ■ x = k± ■ x±_ + k + z = k± ■ x±_ — k + x~ /2, dx — d 2 xj_dz, dk = d z k^dk + and 



We find 



i' a (k) ee J dxe- lk - x e- l ^ A ^u a (x). 



X A (k,P) = (2irf £ -Pa)- 



(2.13) 



(2.14) 



a£F 



The motivation for the 'double-prime' notation is the subject of the next section. 



III. WAVE FUNCTION SUBTLETIES 



It would be useful to express \ A {k, P) m terms of solutions of the ordinary Dirac equation, because one may use a 
standard computer program p8[ . To this end it is convenient rewrite Eq. ( 2.13| ) 



- / dxe- ik " B e- iB ' A ^{x\i/>+). 



(3.1) 



Note that the difference between \ip'a + } of Eq. (3.1) and u" a of Eq. ( 2.13| ) is simply the normalizat ion factor V2P + A. 

These 'double-primed' fields satisfy another versio n of the mode equations Eq. ( |2.6|) and Eq. (2/7) following from 
an application of the Soper-Yan transformation Eq. (2.4), and are given by 



[-id, ~ 9vV+] = [a x ■ (p ± - g v V ± ) + 0(M + gM (3-2) 

[id, + g v V + - 2g v V- + p~] = [a x ■ ip± - 9vV±) + (3(M + g s <j>)] |^"). (3.3) 



If one multiplies Eq. (3.2) by 7 + and Eq. (3^3) by 7 and adds the two equations, using V + = V = V = V , one 
obtains 



[- 7 3 ^ +p-/2) + 7 °(p- 12 - g v V )] ^(x) = [ 7± • p ± + M + g s cf>] ^(x) 



(3.4) 
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Eq. (3.4) is almost the same as the Dirac equation of the equal time formulation (for the ip' fields), with the exception 
of the p~ term multiplying 7 3 . Removing the offending term gives the relationship between the ip' an d ip" fields 



-?„ 2/2 „-ig„A(as) 



x\tp a ) 



(3.5) 
(3.6) 



Eq. ( p.6[ ) is the approximate relationship between the ip and ip' fields determined in Ref. |2l| . The approximation 
lies in the fact that the spectrum condition is not maintained exactly, and the resulting Fourier transform of the 
wavefunction will have unphysical support for k + < 0. This support is largely irrelevant as it manifests far out on an 
expo nen tial tail since p„ co ntai ns the nucleon mass. The relationship between the ip and ip" is exactly our definition 
Eq. (34). The use of Eq. (3.5) in Eq. (3.4) leads immediately to the result that the fields ip' a satisfy the ordinary 
Dirac equation 



7 °(p-/2 - g v V°)iP' a (x) = h-p + M + g s 4>] iP' a (x) 
We now are ready to derive a representation of Eq. ( |2.l| ) in terms of these nucleon wave functions. 



(3.7) 



IV. DERIVATION OF THE PLUS MOMENTUM DISTRIBUTION 

In Ref. , it was determined that a plus momentum distribution in QHD-I is given by 

/(*+) = 2 £ I d 2 x ± \(x ± ,k+\^)\ 2 ■ 



(4.1) 



This distribution peaks at k + /M = y ~ 0.8 for 16 0, (with smaller values for heavier nuclei) but is not the distribution 
obtained from the covariant formalism of Section II. The connection between this f(y) and the covariant /at(?/) was 
made in Ref. [ fl7j |; it was determined that, in the limit of infinite nuclear matter, the relationship between f(y) and 
/n(jj) is simply a shift in the argument by the vector meson potential: 



f(y) = f N (y + 9vV + /M). 



(4.2) 



This shift arises from the use of the Soper-Yan transformation Eq. (2.4) where the ip' fields are those appearing in the 
Lagrangian and are used to determine /j\r(j/), whereas the ip fields are used to determine f{y). In finite nuclei, this 
relationship is som ewha t more complicated since the vector meson potential is no longer a constant over all space. 
We start with Eq. ( 2.14 ), and see that 



Tr~,+ X A (k,P) = (2^) 2 ^Tr[ 7 +<(feK(fc)]^(fc--p-) 
= 16^P+ J 4]T|(fcK+)| 2 : S(k- -p~). 



Substituting into Eq. (2.1) we obtain 

f N (v) 



J2 I dk5(y-k+/M)\(k\^)\ 

aeF 



Use of Parseval's identity and integrating over k + gives us our main result: 



d 2 x ± \(x ± ,My\^+) 



(4.3) 



(4.4) 



so the p lus momentum distribution is related to Fou rier transform of the ip" wave functions. One can see the similarity 
to Eq. (4.1); the difference lies entirely in Eq. (3.1). It should be emphasized that this result does not depend on the 
approximation in Section III. 

We shall use Jn{v) to compute the nuclear structure function F2a(%) in Section |v|, b ut first we derive a version of 
the Hugenholtz-van Hove theorem valid for finite nuclei. To do that, multiply Eq. (4.3) by y and integrate 



(y) = I dyyf N (y) 



(4.5) 



a£F ' 
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Now remove the plus projections and re-express if>" and its complex conjugate in coordinate spaces x and x' . One 
can then integrate over k yielding a delta function S(x — x') which allows integration over x' 



(y) 



= £ / ^< t (^)7°7 + ^ + C(= 



We wish to look at the ip fields in order to understand our result in the context of Ref . |^T| , so we need to perform 
the Soper-Yan transformation Eq. and use X^j° = X 



(y) 



= E 

M 



dx^ a (x)-y + [id + + g v V + (x)]ip a (x) 



If we explicitly put in the the nuclear state vectors, we can perform the sum on a by inserting creation and annihilation 
operators; we can add the time dependence for free since it is unaffected by d + an d ca ncels with both fermion fields, 
and the vector potential is static. We have effectively undone the substitution Eq. ( |2.5| ) and now have an expectation 
value of an operator 



(4.6) 



Using the vector meson field equation in QHD-I 

8^+ + mlV + = gjrf + il>, 



integrating by parts, and anti-symmctrizing one can re-express the second term of Eq. (4.6) 



(y) 



J dx(i;j + id + ijj + m 2 v V + V+ + V +fl V fl + ) 
dx(T ++ - d + (/)d + (t>} 



1 

Ma 

— (P+-P+) 
Ma [ s ' 

P+ 

M A 



(4.7) 



where T ++ is the canonical energy momentum tensor, is the plus momentum of the scalar meson fields, and P + 
is the total nuclear plus momentum. The result Eq. (4.7) constitutes an analog of the Hugcnholtz-van Hove theorem 
p9| for finite systems; the equality becomes exact in the nuclear matter limit, where the scalar meson contribution 
vanishes, as shown in our previous work j^j]. This means that we may anticipate that the binding effect will again be 
small. The 'mixin g' of the vector operators and the scalar meson contribution will be elaborated on in a more general 



context in Section VI 



It is also worthwhile to explicitly evaluate the expression Eq. (4.4) for /^(y) in the limit of infinite nuclear matter. 
In this case, V° = V + = V~ are constant and V± = 0, so we find 



A(z, x 



dz'V°(z',x ± ) 



-V°z 
-V+z, 



(4.8) 



so that Eq. (3.1) becomes 



dxe 



-i(k+-g v V + )z i 



Therefore Eq. (4.4) becomes 



f N (y) 



2M^2 ld 2 x x \(x x ,My-g v V+\TP+)\ 2 , 



(4.9) 



(4.10) 



which is simply the expression (4.1) modified by a shift in the argument of g v V + / M. Thus we find Eq. (4.2) is satisfied 
in the nuclear matter limit. It is important to stress that all that is recovered here is the shift in the argument and 
not any particular form of the plus momentum distribution which arises from the specific model used. 
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FIG. 2: Plus momentum distributions, /jv(y), for 40 Ca (solid), le O (short dashes), 208 Pb (dot-dashes) and nuclear matter (long 
dashes) . 




FIG. 3: f(y) distributions for 16 and nuclear matter after application of the Soper-Yan transformation. Note that the peaks 
occur at y < 1. 



V. NUCLEAR STRUCTURE FUNCTIONS 



We determine the wave functions appearing in Eq. (4.4) numerically from a relativistic self-consistent treatment 
following Horowitz and Serot fl27| using the same program j28| which includes electromagnetic effects. The plus 
momentum distribution follows and is given in Fig. g for 16 0, Ca, 208 Pb and in the nuclear matter limit (the 16 
calculation is also s how n in Fig. ^) . One can see that the peaks appear near y = 1 as required by the Hugenholtz-van 
Hove theorem Eq. (4.7). 



It is worth noting that application of the Soper-Yan transformation Eq. (2.4) to the ip" wavefunctions obtained 



from the equal time wavefunctions reproduces the plus momentum distributions, inclu ding the correct asymmetry, 
of the light front calculations in Ref. ^J, which did not use the approximation Eq. (3.5), as shown in Fig. || for 
Oxygen and nuclear matter. This demonstrates the excellence of the approximation relating light front and equal 
time wavefunctions. One can see that the effect in finite nuclei of the Soper-Yan transformation is to shift and broaden 
the plus momentum distribution, w hile in nuclear matter it is just a shift. If these distributions wer e to be used in 
the nuclear structure function Eq. (5.2) though, since (y) ~ 0.8 for Oxygen, the ratio function Eq. (5.4) would fall 
precipitously to nearly zero at x ~ 0.6 in stark contradiction with experiment. 



We also evaluated the plus momentum distribution Eq. (4.4) with the simple non-relativistic harmonic oscillator 
shell model as an additional check on our method. These (equal time) wavefunctions give us an explicit, although 
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FIG. 4: Plus momentum distribution for 16 O calculated with harmonic oscillator (solid curve) and relativistic Hartree (dashed 
curve) wave functions. 



approximate in the sense of Section III, closed form of the plus momentum distribution for b O 



f N (y) ~ t[ e -?<n.-v) a 

V l07T 

+2(i + e(m P -y) 2 )e~^- y ~>\ (5.1) 

with r)ni = M + w(2(n - 1) + I + 3/2) - v where v ^ 50 MeV, and £ = Jib where 6 = (mi)" 1 ' 2 ~ 1.6 fm is the 
oscillator length which is fit to the root mean square radius of Oxygen (i? 2 ) 1 / 2 ~ 2.7 fm. The distribution Eq. (pjl) 
narrows for larger £ which corresponds to an increasing root mean square radius. This distribution is plotted in Fig. || 
where one can see that it peaks near y = 1 like the relativistic Hartree calculation, but appears to have a smaller 
value of (y) . It is worth noting that the Hartree calculations are in the relativistic equal time framework and put into 
our relativistic light front formalism, while the harmonic oscillator calculations are non-relativistic and put into our 
relativistic formalism. 

The structure function is given by the convolution 

= £ dyf N (y)F 2N (x A /y), (5.2) 

with xa = Q A/2 P ■ q = xM/M. The assumption that nuclear effects do not modify the structure of the nucleon is 
embodied in Eq. ( |5.2| ) by the use of the structure function of a free nucleon; we use the parameterization 

F 2N {x) = 0.58^(1 - x) 2 - 8 + 0.33V^(1 - xf* + 0.49(1 - x) 8 . (5.3) 

The experiments measure the ratio function, defined as 

*(*) = ^fy (5-4) 

The results of our calculations are plotted for 16 0, 40 Ca, 208 Pb and in the nuclear matter limit in Fig. |5| showing data 
for Carbon, Calcium and Gold from SLAC-E139 |[o) and an extrapolation |3^| for the nuclear matter calculation. 
The most striking result is that these calculations fail to reproduce the EMC effect; the curves consistently miss the 
minima in the data, and the agreement gets worse with increasing A. Another important result is that the ratio 
function does not fall to zero as would be the case if the small effective mass (~ 0.56A/ for n ucle ar matter in QHD-I) 



were the relevant parameter describing the binding effect which would follow from using Eq. (4.1) instead of Eq. ( fh4| ). 
The results also show a minimum near x ~ 0.6 for Oxygen and nuclear matter that is deeper than the Calcium and 
Lead calculations. This is a curious feature that contradicts the trend in experimental data, and is due to the effects 
of two parameters. 

The first, and most important, is that of the location of the peak of the plus momentum distribution given by 



Eq. (4.7), which gradually approaches y = 1 as the nuclear matter limit is reached. This is due to the fact that scalar 
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FIG. 5: Ratio functions for 16 0, 40 Ca and 208 Pb showing data for Carbon, Calcium and Gold, respectively, from SLAC-E139 
[BOl. The nuclear matter calculation shows extrapolated data [B2[ . 



mesons carry a small amount of plus momentum |21| that vanishes as A — > oo. The closer to y — 1 the peak is in 
Fig. the less pronounced the minimum in Fig. |^, all else remaining constant. The second effect is due to M, which 
reaches a minimum at 56 Fe corresponding to a more pronounced minimum of the ratio function than for A < 56 or 
A > 56, keeping the scalar meson contribution constant. 

Using a Gaussian parameterization of the plus momentum distribution and the experimental binding energy per 
nucleon via the semi-empirical mass formula, we have modeled the dependence of the minimum of the ratio function, 
R(x ~ 0.72), on the number of nucleons in the nucleus in Fig. ^. The motivation for the use of Gaussian plus 
momentum distributions is based on the expansion @] 



F 2 a(xa) = F 2N (xa) + tx A F2 N 
+-l[2x A F! 2N (x A ) 

where e = 1- J dyyf N {y) 



-r 2 F" 
A 2N 
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dy(y 



i) 2 f N {y) 



(5.5) 
(5.6) 

(5.7) 



The Gaussian parameterization uses the peak location and width, (y) and ((y 2 ) — (y) 2 ) 1 ^ 2 respectively, from the 
relativistic Hartree calculations in Fig. ^|, and is normalized to unity. This allows us to obtain a plus momentum 
distribution for any A with minimal effort. We show the combined effect of scalar mesons and binding energy per 
nucleon on the ratio function along with the effect of scalar mesons alone using a constant binding energy per nucleon 
of —8.5 MeV independent of A. It can be seen that a changing M with A has the most effect for nuclei much larger 
than Iron, but does not change the general trend that the minimum of the ratio function becomes less pronounced as A 
increases due to the vanishing scalar meson contribution and the peak of the plus momentum distribution approaching 
unity. This dependence of the binding effect on A is quite different, both in magnitude and shape, than the trend 
in experimental data summarized in Ref. p^] which satisfies R{x ~ 0.72) ~ A -1 / 3 , so that the minimum becomes 
more pronounced as A increases. This fully demonstrates the inadequacy of conventional nucleon-meson dynamics to 
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FIG. 6: R(x — 0.72) as a function of A including scalar meson and binding effects (solid line), and leaving binding energy per 
nucleon constant at —8.5 MeV (dashed line). The data are from SLAC-E139 |30|. 



explain the EMC effect. 



VI. SCALAR MESON CONTRIBUTION TO PLUS MOMENTUM AND MORE GENERAL 

CONSIDERATIONS 



The average value of y, given by Eq. (4/7), yields the nucleon contribution to the plus momentum, and is less than 
one which can be seen in Fig. | We now address the remaining plus momentum in finite nuclei. Previous results 
fHf show that a small fraction (5y ~ 0.005) of the plus momentum is carried by the scalar mesons which vanishes as 
the nuclear matter limit is approached. This is due to the fact that scalar mesons couple to gradients in the scalar 
density (arising mainly from the surface of finite nuclei) which vanish as A — > oo. The question is: why are scalar 
mesons allowed to carry plus momentum and not vector mesons? 



The simplest answer lies in the Dirac structure of Eq. (2.1); the 7 + in the trace picks out terms in the full 
interacting Green's function with an odd number of gamma matrices which includes all Lorentz vector interactions 
and excludes Lorentz scalar interactions. The Dirac structure of /jv(y) is directly related to the Dirac structure of the 
energy momentum tensor, so the answer also lies there and illuminates a problem with conventional nucleon-meson 
dynamics. The component of the energy momentum tensor relevant to the plus momentum, from a chiral Lagrangian 
containing isoscalar vector mesons, scalar mesons and pions, is given by Jl8|, [l9| 

T++ = V + ^V^ + + m 2 v V + V + +Jn + id + ij) + d+cj)d + (j) 

+d+n ■d+n + Tv-d+n 7 '' d, +?r (I-/! S in 2 (6.1) 

TX A 7T J 

Since each of the terms in Eq. ( |6.1[ ) involves one of the fields, it is natural to associate each term with a particular 
contribution to the plus momentum. This decomposition, though, is not well defined; field equations relate various 
components. We see the first three terms of Eq. ( |6.1[ ) appear in (y), which defines the nucle on contribution to the 
total nuclear plus momentum, in the derivation of the Hugenholtz-van Hove theorem Eq. ( |4.7[ ) ; we are not allowed to 
have the vector mesons contribute a well defined fraction of plus momentum. This means that one could trade certain 
mesonic degrees of freedom for nucleons by replacing mesonic vertices with nucleon point couplings, for example, in 
line with the general concept of effective field theory. In our case the first three terms are related by the vector meson 
field eq uati on, but the fourth is left out since the scalar mesons couple to the scalar density ipip which is not present 
in Eq. ([0]). Therefore the scalar mesons (and pions) contribute a well defined fraction of plus momentum. These 
explicit meson contributions create an EMC binding effect, but the pionic contributions are also limited by nuclear 
Drell-Yan experiments |36]| to carrying about 2% of the plus momentum which is insufficient to account for the entire 
EMC effect which corresponds to about 5% of the plus momentum for Iron. 
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VII. SUMMARY AND DISCUSSION 



The minimum in the EMC effect is known to have a monotonically decreasing behavior with A, which has been 
studied in Refs. ||(], Q among others. Our present theory is defined by the use of the mean-field approximation, 
along with the assumption that nuclear effects do not modify the structure of the nucleon. This theory leads to 
results in severe disagreement with experiment. Not only do we find that the depth of the minimum is monotonically 
decreasing with A, but it has a smaller magnitude than experiment. These results, which fail to capture any of the 
important features of the experiments, represent a failure of relativistic mean field theory. Furthermore, the plus 
momentum distributions we compute give (y) ~ 1 which indicates that nearly all of the plus momentum is carried by 
the nucleons. In order to reproduce the data, the nucleon plus momentum must be decreased by some mechanism that 
becomes more important at larger A. Nucleon-nucleon correlations cannot take plus momentum from nucleons, and 
explicit mesonic components in the nuclear Fock state wavefunction carrying plus momentum are limited 33 , 34|. |3^] 
by Drell-Yan experiments |3(| . Thus it appears that the EMC effect may be due to something outside of conventional 
nucleon-meson dynamics. For example, true modifications to nucleon structure caused by nuclear interactions could 
be important, in which case one would need to use models such as the mini-delocalization model quark- meson 
coupling (QMC) model p?], 38, 39 or the chiral quark soliton model reviewed in Eof to include those effects. 
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